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ABSTRACT 

We present a measurement of the angular power spectrum of the cosmic microwave background 
(CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution 
made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 deg 2 
of the southern sky, in a 4? 2- wide strip centered on declination 53° South. The CMB at arcminute 
angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel’dovich (SZ) effect 
from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 
brightest point sources in our maps, we estimate the power spectrum between 600 < £ < 8000 using 
the adaptive multi- taper method to minimize spectral leakage and maximize use of the full data set. 

Our absolute calibration is based on observations of Uranus. To verify the calibration arid test the 
fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and 
recover the WMAP power spectrum from 250 < £ < 1150. The power beyond the Silk damping tail 
of the CMB (£ ~ 5000) is consistent with models of the emission from point sources. We quantify the 
contribution of SZ clusters to the power spectrum by fitting to a model normalized to as = 0.8. We 
constrain the model’s amplitude A$ % < 1.63 (95% CL). If interpreted as a measurement of as, this 
implies af 2 < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to 
a 6-parameter ACDM model plus point sources and the SZ effect is consistent with these results. 

Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

The cosmic microwave background (CMB) radiation 
captures a view of the universe at only ~ 400, 000 years 
after the Big Bang. The angular power spectrum of 
temperature anisotropies in the CMB has been crucial 
in developing the current standard cosmological model, 
in which the universe today contains some 5% baryonic 
matter, 23% dark matter and 72% dark energy. We re- 
fer to this model throughout as the ACDM model. The 
temperature power spectrum has been measured to good 
precision for multipole moments £ < 3000. At angular 
scales with £ < 2000, the power spectrum matches the 
predictions of ACDM, and it can be used to constrain 
multiple parameters of the cosmological model (e.g., 
Dunkley et al. 2009; Brown et al. 2009; Reichardt et al. 
2009b; Sievers et al. 2009). The agreement between 
the current polarization anisotropy measurements and 
the predictions of ACDM cosmology further supports 
the model (e.g., Dunkley et al. 2009; Brown et al. 2009; 
Chiang et al. 2009). 

At £ > 3000, the signal from the primary CMB 
anisotropy becomes dominated mainly by two popu- 
lations. The first consists of point source emission 
from both radio and dusty infrared-emitting galaxies. 
The second is the population of massive galaxy clus- 
ters that give rise to the Sunyaev-Zel’dovich (SZ) ef- 
fect (Sunyaev & Zel’dovich 1970), in which CMB pho- 
tons scatter off the electrons of the hot intra-cluster 
medium. Removal of the brightest foreground galax- 
ies and SZ clusters can reduce their contributions to 
the power spectrum, but the net anisotropy power due 
to unidentified sources still dominates the exponentially 
falling primary anisotropy spectrum at small scales. 

The top panel of Figure 1 shows recent measurements 
of the anisotropy at £ > 2000. In analyses of 30 GHz 
data, radio point sources are masked out and a resid- 
ual component is modeled and subtracted. For 150 GHz 
analyses, radio sources are masked and the residual radio 
and dusty galaxy contribution is estimated. 

The SZ effect has a unique frequency signature. In 
CMB temperature units (the units of Figure 1) the am- 
plitude at 150 GHz is roughly half that at 30 GHz. The 
amount of SZ power is governed by erg, which measures 
the amplitude of the cosmic matter power spectrum on 
8h“ 1 Mpc scales; the SZ power scales approximately as 
<rj(0&h) 2 (Seljak et al. 2001; Komatsu Seljak 2002). 
At 30 GHz, CBI reports (Sievers et al. 2009) excess emis- 
sion above the ACDM model at £ & 3000 after account- 
ing for all known radio sources (Mason et al. 2009). A 
possible source of the “CBI excess” is the SZ effect. 
The SZA data (Sharp et al. 2009), also at 30 GHz and 
also after accounting for point sources, are consistent 
with the ACDM model at £ ~ 4000. The ACBAR re- 
sults (Reichardt et al. 2009b) at 150 GHz are consistent 
with the CBI excess and ACDM. Sievers et al. (2009) 
show that all the above data in the top panel of Figure 1 
are consistent, within 95% CL, with ACDM plus a SZ 
contribution of erg = 0.922 =t 0.047 (la error bars). Re- 
cently, the South Pole Telescope (SPT) group reported 
<rg == 0.773 ± 0.025 based on the power spectrum at 150 
and 220 GHz (Lueker et al. 2009). 

In this paper, we present a new measurement of the 
CMB anisotropy power spectrum in the range 600 < £ < 


8000 (corresponding to angular scales of approximately 
1.4' to 18'). The observations were made at 148 GHz in 
2008 with the Atacama Cosmology Telescope. Figure 1 
(bottom panel) shows the results, which are discussed in 
Section 7. We briefly describe the instrument, the obser- 
vations, the calibration of the data, and the map-making 
procedure. The dynamic range of the power spectrum 
is large enough that we employ new techniques in spec- 
tral estimation. We confirm that difference maps of the 
data give power spectra consistent with no signal and 
that the power spectrum is insensitive to several details 
of our analysis. We use the shape of the power spectrum 
to bound the dusty galaxy contribution and the power 
from SZ clusters. 

2. INSTRUMENT AND OBSERVATIONS 

The Atacama Cosmology Telescope (ACT) is 
a 6-meter, off-axis Gregorian telescope optimized 
for arcminute-scale CMB anisotropy measurements 
(Fowler et al. 2007; Hincks et al. 2009). It was installed 
at an elevation of 5190m on Cerro Toco 1 in the Atacama 
Desert of northern Chile in March 2007. Observing 
conditions in the Atacama are excellent owing to the 
elevation, the arid climate, and the stability of the at- 
mosphere. After all cuts, the median precipitable water 
vapor (PWV) was 0.49 mm during the observations 
presented here. 

The Millimeter Bolometer Array Camera (MB AC), 
the current focal-plane instrument for ACT, uses high- 
purity silicon lenses to reimage sections of the Grego- 
rian focal plane onto three rectangular arrays of detec- 
tors. The arrays each contain 1000 transition edge sensor 
(TES) bolometers. Their spectral coverage is determined 
by metal- mesh filters (Ade et al. 2006) having measured 
band centers of 148 GHz, 218 GHz, and 277 GHz. The 
bolometers are cooled to 300rnK by a two-stage helium 
sorption fridge backed by commercial pulse-tube cry- 
ocoolers. The telescope performance and control sys- 
tems (Hincks et al. 2008; Switzer et al. 2008), camera de- 
sign (Thornton et al. 2008; Swetz et al. 2008), detector 
properties (Zhao et al. 2008; Niemack et al, 2008), and 
readout electronics (Battistelli et ah 2008) are described 
elsewhere. 2 

2.1. Observations 

ACT operated with all three arrays from mid- August 
to late December, 2008 for the data presented here. The 
observing time was divided between two regions away 
from the galactic plane. The deepest observations cover 
900 deg 2 of the southern sky in a strip 8° wide centered 
on 6 = -53° with RA from 19 h to 24 h and 0 h to 7 h 36 m . 
The analysis presented in this paper uses only data from 
the central 228 deg 2 of the southern strip and only ob- 
servations made with the 148 GHz array. The slightly el- 
liptical beam has full-widths at half- maximum (FWHM) 
of 1.40' by 1.34' at this frequency (Hincks et al. 2009). 

The observations were made by scanning the sky at a 
constant elevation of 50°. Each scan is 4? 5 wide on the 
sky (7?0 in azimuth angle). Scans repeat every 10.2s. 
Each half-scan consists of 4.2 s of motion at a constant 

1 ACT is at 22.9586° south latitude, 67.7875° west longitude. 

2 The site http://www.physics.princeton.edu/act/ archives 
papers by the ACT collaboration. 



4 


J. Fowler et al. 



20 200 500 1000 1500 2000 3000 4000 5000 6000 7000 8000 



20 200 500 1000 1500 2000 3000 4000 5000 6000 7000 8000 

Multipole l 


Fig. 1. — Recent measurements of the CMB power spectrum, including this work. Top : the measurements of WMAP (Nolta et al. 
2009), Bolocam (Sayers et al. 2009), QUaD (Brown et al. 2009; Friedman et al. 2009), APEX-SZ (Reichardt et al. 2009a), ACBAR 
(Reichardt et al. 2009b), SZA (Sharp et al. 2009), RIMA (Dawson et al. 2006), CBI (Sievers et al. 2009), and SPT (Lueker et al. 2009). For 
all the results, a radio point source contribution has been removed either by masking before computing the power spectrum (at 150 GHz), 
or by masking and modeling the residual (at 30 GHz and for WMAP). APEX-SZ additionally masks clusters and potential IR sources. 
Bottom: The ACT power spectrum from this work. The inset shows the cross-power spectrum between ACT and WMAP maps in the AGT 
southern field (see Section 3.5), which we use to check both the validity of the maps at larger scales and the absolute calibration. Only the 
ACT power spectrum is analyzed in this paper. In both panels and the inset, the solid curve (blue) is the ACDM model of Dunkley et al. 
(2009) (including lensing). The SZ effect and foreground sources are expected to contribute additional power, as shown in Figure 4 and 
Table 1. For display purposes- — and only in this figure — we scale our result by 0.96 in temperature relative to the Uranus calibration; this 
calibration factor best fits our data to the ACDM model and differs from the Uranus calibration by O.Tcr. Recent WMAP observations of 
Uranus suggest the same rescaling factor (see footnote to Section 3.2). ACT bandpowers for t > 4200 have been combined into bins of 
A£ = 600 for this figure: they are given in a note to Table 1. 

speed of l?5/s followed by 0.9 s of acceleration. The first apart from small areas at the edges, 
half of each night is spent observing the field rising in the As the telescope scans in azimuth at constant elevation, 

eastern sky, after which ACT turns to the western sky each detector is sampled at 399 Hz. The data sampling, 
to observe the same field as it sets through the standard position reading, and all housekeeping data are synchro- 
elevation of 50° . The scan strategy is designed to mini- nized by a shared 50 MHz clock; absolute times are ref- 

mize changes in the telescope's orientation with respect erenced to a GPS receiver with 0.25 ms accuracy. The 

to the local environment while ensuring cross-linked oh- data are stored in continuous fifteen-minute segments 

servations in celestial coordinates. Sky rotation ensures called time-ordered data sets (TODs). Each TOD re- 

that all detectors sample all points in the field each night, quires 1.6 GB of storage per detector array, or 600 MB 
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after applying lossless compression. 

In addition to the main survey, we perform occasional 
calibration measurements. Most nights, a few minutes 
are used to measure a planet when it passes through our 
standard elevation. Any given target is used only once 
in every three nights to minimize non- uniformity in the 
coverage of the CMB regions. The planetary observa- 
tions allow us to measure the system’s relative and abso- 
lute responsivities, pointing, and beam profiles; they also 
provide a way to check for any time variations in each. A 
series of tuning and biasing procedures are also followed 
each night before and after regular CMB observations, 
to optimize and roughly calibrate the detector response. 

3. DATA REDUCTION AND MAP-MAKING 

The goal of the data reduction is to estimate the 
maximum- likelihood map of the sky. We select properly 
tuned detectors and calibrate their pointing and relative 
gains. The absolute gains are determined by observa- 
tions of Uranus. The map-making itself is an iterative 
and computationally intensive process. 

3.1. Data selection 

The first step in the data reduction is to select TODs 
when the receiver is cold enough for stable operation and 
when the precipitable water vapor (PWV) is less than 
3 mm. Of the « 2880 h between Aug 25 and Dec 24, 
2008, MBAC was on line for 1352 h, roughly correspond- 
ing to the night time hours. After cutting on PWV and 
instrument performance, 1031 h (35% of calendar time) 
went into the pipeline. Of those, 850 h were spent ob- 
serving the southern strip. 

In the second step of the data reduction, we check 
whether each detector’s data exhibit problems that war- 
rant removal from the final analysis (Dixnner 2009). In 
the case of rare and transient effects, data are removed 
for a single detector up to a few seconds at a time; other 
problems are diagnosed in 15-minute intervals, and when 
necessary, each detector’s data may be removed for the 
rmtire TOD. We check the detector’s response speed; look 
for transient effects such as cosmic ray hits; and en- 
sure that the feedback loop remained locked, keeping the 
SQUID amplifiers in the linear regime (Battistelli et al. 
2008). Finally, we compare the detector’s response to 
atmospheric emission with that of the array- wide aver- 
age, in each case filtering out data at frequencies above 
50 mHz. Detectors are cut if their data do not have at 
least a 0.98 correlation with the average or if the ampli- 
tude of the atmospheric signal is not within 15% of the 
median amplitude found among good detectors. This 
last cut is the largest and typically removes 100 detec- 
tors. Additionally, we bin the data by scan and compute 
the rms, skewness, and kurtosis in each bin. A cut on 
rms removes an additional 30 detectors on average, while 
the combined cut on skewness and kurtosis removes only 
one additional detector. When the instrument is well- 
tuned, ACT typically has 680 science- grade detectors at 
150 GHz. The remainder are discarded from the analy- 
sis or used to monitor instrumental effects. The effective 
sensitivity of the array, including all instrumental and at- 
mospheric effects, is ~ 30 pK s 1//2 in CMB temperature 
units. 

3.2. Calibration to planets 


After we have a reliable set of detectors to use, the 
next step is gain calibration. In practice, the cuts are 
computed twice, so that the final set of cuts is based on 
calibrated data. We use Uranus as our absolute reference 
standard. 

The atmosphere is an excellent continuous flat-fielding 
calibrator. We find that the time variation in detector 
relative responsivities is not significant. Therefore, the 
detectors’ relative responses to the atmosphere (averaged 
over the entire season) are used to convert each output to 
a common scale. Observations of Saturn confirm these 
relative calibrations; it is the only source observed in 
2008 bright enough to use for this confirmation. 

We take the brightness temperature of Uranus at 
148 GHz to be Tu — 112±6K (Griffin &, Orton 1993; 
Marten et al. 2005; Kramer et al. 2008). All results in 
this work scale in proportion to this temperature. For 
planet calibrations, we use the beam solid angle = 
218 ± 4nsr (Hincks et al. 2009). The net calibration er- 
ror of 6% in temperature is the combination of these 
two uncertainties and is dominated by uncertainty in the 
brightness of Uranus. 

Saturn was used as a rough check on the absolute cal- 
ibration result. The Saturn observations suggest that 
the ACT data are 10% brighter than inferred from the 
Uranus data. We attribute this inconsistency to the dif- 
ficulty in modeling the brightness of the planet and its 
rings as they vary over time and across frequencies. 

We track the stability of the system over time in a num- 
ber of ways. The absolute celestial calibration is checked 
through measurements of a planet on most nights. The 
conversion between raw data and units of power absorbed 
on the detectors is calibrated in two ways. The conver- 
sion is estimated twice nightly by modulating the TES 
detector bias voltage with a small additional square-wave 
and measuring the response. A calibration is also per- 
formed at the beginning of each night by sweeping all 
detectors through the full range of bias voltages. Both 
methods are described further in Niemack (2008) and 
Fisher (2009). Based on several tests (Switzer 2008), we 
find that the gain of each detector is constant over the 
season to better than 2%. 

As part of the calibration, the detector temporal re- 
sponse is deconvolved from the detector time streams. 
We model the response with a single time constant. The 
detectors’ time constants are determined from the planet 
measurements and bias steps (Hincks 2009). The median 
time constant is r c i e t. = 1.9ms (thus / 3 C ib = 84 Hz), which 
corresponds at the ACT scan rate to t ~ 31,000. Slow ? 
detectors with f'sdB <15 Hz are not used. We also decon- 
volve the anti-aliasing filter imposed in the data acqui- 
sition system (Battistelli et al. 2008). We then further 
filter and sample the time stream at half the raw rate to 
speed the map-making step. This last anti-aliasing filter 
acts only at angular scales smaller than the beam and 
does not affect the maps. 

3.3. Pointing reconstruction 

The planet Saturn is bright enough that each bolome- 
ter detects it with high signal- to-noise in a single 
scan. We fit in the time domain to find the best two- 
dimensional location for each detector relative to the 
notional array center. The fit produces altitude and 
azimuth offsets for each, along with information about 
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Fig. 2. — The map and a difference map of the ACT southern field at 150 GHz. The same filters used in the power spectrum analysis 
are applied to both maps: an isotropic high-pass filter suppresses power for £ < 300 (Equation 2), and all modes with \£ x \ < 270 are set to 
zero, as described in Section 5.2. Top : The ACT southern field. The intensity scale is pK (CMB units). The 360 deg 2 with lowest noise 
are shown. The squares 4? 2 on a side indicate the 13 patches used for the CMB power spectrum analysis (228 deg 2 total). Middle : A 
difference map made from two halves of the same data set. Most of the remaining structure visible at large scales is well below the range 
of £ that we consider in the power spectrum analysis. Bottom: The rms temperature uncertainty for one square arcminute pixels. 


the relative gains, the beam sizes, and the detector 
time constants. The detector pointings are consistent 
with optical models of the telescope and reimaging op- 
tics (Fowler et al. 2007). The relative pointing used in 
this work is an average over 22 observations of Saturn. 
The rms uncertainty is 1.2" for the detector relative 
pointings, which contributes negligibly to any pointing 
error in the sky maps. 

The location of the array center can also be found from 
each observation of Saturn. It has a scatter of 4.3" rms 
over the 22 observations, which we attribute to slight 
thermal deformations in the telescope mirrors and their 
support structures. 

The Saturn data give the absolute location of the ar- 
ray center only for other observations taken at the same 
horizon coordinates. To determine the absolute location 
of the array center during science observations, we used 
approximately 20 known radio sources found in prelimi- 
nary maps. We find this approach both simpler and bet- 
ter constrained than making a complete pointing model 
of the telescope. The ACT southern field was observed 
both rising and setting, at azimuths centered at 30° on 
either side of south. These preliminary maps were made 
separately for rising and setting data and were used to 
determine separate pointing corrections for the two cases. 
After correction, we estimate the maps to have 5" point- 
ing uncertainty. 

3.4. Map - making 


The goal of the map-making step is to take the 3200 GB 
of cut, calibrated, and deconvolved raw data from both 
the rising and setting scans and produce from it a maxi- 
mum likelihood estimate of the sky. This work produces 
a 200 MB map, 16,000 times smaller than the raw data 
set. We first multiply the data for each detector in each 
TOD by a window function, reducing the weight in the 
first and last 10 seconds of each file, then remove a single 
offset and slope for the entire 15-minute period. There 
is no additional filtering in the time domain, though the 
lowest frequencies are given no weight in the process of 
maximizing the likelihood of the final map. The detector 
data are combined into a data vector d. The mapping is 
done in a cylindrical equal-area projection with a stan- 
dard latitude of S = —53.5° and pixels of 30" x 30", 
roughly one third of the beam FWHM. The map is rep- 
resented as a vector rn of length jY p * x ~ 10 r . We model 
the data as d = Pm + n where the matrix P projects 
the map into the time stream, and n is the noise, which 
has covariance matrix N. The maximum likelihood solu- 
tion, m, is given by solving the mapping equation (e.g., 
Tegmark 1997): 

P T N“ 1 Pm - P T N -1 dL (1) 

We solve for fri iteratively using a preconditioned 
conjugate gradient (PCG) method (Press et al. 2007; 
Hinshaw et al. 2007). Based on simulations, we find 
that the solution is an unbiased estimator of the sky for 
£ > 600, the multipoles that we analyze iri this paper. 
Additionally, the clear cross-correlation with WMAP 
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(Section 3.5) indicates that the maps are likely unbi- 
ased as low as £ « 200. This approach is different 
from that taken by other recent measurements of the 
fine angular scale anisotropy in which the time stream 
data are filtered and binned (e.g., Friedman et ah 2009; 
Reichardt et al. 2009a; Lueker et al. 2009). 

Solving Equation 1 requires careful consideration of the 
noise structure of the ACT data. First, the data are 
weighted in Fourier space to account for the variation in 
noise with frequency, particularly to reduce the noise that 
the atmosphere adds at low frequencies. For frequencies 
less than 0.25 Hz, we set the statistical weights to zero. 
Higher frequencies are given successively more weight in 
inverse proportion to the noise variance in each band. 
Second, we find that there are several modes, particular 
combinations of the 1000 detectors, that correspond to 
signals unrelated to the celestial temperature; for exam- 
ple the common mode, atmospheric gradients, and other 
detector correlations. We handle them by finding the 
ten modes with the largest eigenvalues in each TOD and 
solving for their amplitudes as a function of time along 
with the map m. The power spectrum (Section 6) is 
robust to halving or doubling the number of modes. 

Maps are made using only fractions of the valid data 
so that null tests can be performed. The subset maps 
are also used for finding the CMB power spectrum, as 
described further in Section 5.2. The final map used for 
source-finding, however, is a complete run based on all 
the valid data. Figure 2 shows the map and the region 
used. The computational task is considerable. One iter- 
ation takes 100 seconds on 5000 cores of Canada’s SciNet 
GPC cluster each running at 2.53 GHz; computations re- 
quired before the first iteration take an additional time 
equivalent to approximately twenty iterations. A con- 
verged map requires hundreds of iterations and approxi- 
mately ten CPU years. 

3.5. Calibration to WMAP 

The absolute calibration described in Section 3.2 is 
limited by the uncertainty in the 148 GHz brightness 
of Uranus, the primary calibrator. We complement the 
planetary calibration by cross-correlating with the W- 
band (94 GHz) measurements of WMAP (Hinshaw et al. 
2009). In the angular scale range of 200 < £ < 1200, 
both instruments measure the sky with sufficient signal- 
to-rxoise to permit the comparison. 

The technique requires a single map from WMAP and 
two maps with independent noise from ACT. We assume 
that only one relative calibration ratio between ACT and 
WMAP must be estimated. The WMAP data set con- 
sists of the high-resolution 3 W-band five-year map. The 
two independent ACT maps are each made with one half 
the current data set. We compare the maps only in the 
228 deg 2 strip used for the present power spectrum anal- 
ysis (described in Section 5.1), where the noise is low and 
uniform. 

We compute those cross- spectra in two-dimensional £- 
space that combine either the full ACT data set with 
the WMAP map or the two ACT maps with each other. 
We average over the polar angle in t-space to get a one- 
dimensional cross-spectrum for each. In the case of the 

3 The WMAP maps are at HEALPix resolution N s - lc { e — 1024, 
with 3.5' pixels. 


ACT spectra, the noise is not isotropic, so the spectra 
are angle- averaged with appropriate weights. We then 
find the single calibration factor that minimizes y 2 for 
a model in which the minimum- variance weighted com- 
bination of the ACT- ACT and the ACT- WMAP cross- 
spectra equals the WMAP all-sky spectrum. The model 
accounts for the fact that the ACT and WMAP mea- 
surements have noise varying as different functions of £. 

The result is a calibration factor with < 6% fractional 
uncertainty in temperature. The calibration factors de- 
rived from WMAP and from Uranus observations are 
consistent to < 6%. 

4. FOREGROUNDS 

Radio and infrared galaxies are the dominant sources 
of foreground emission in the 150 GHz band at £ > 
1000. To study the underlying CMB spectrum, the 
sources must be identified and masked, and resid- 
ual contamination in the power spectrum must be 
accounted for. Many approaches to the problem 
have been described in the CMB literature (e.g., 
Wright et al. 2009; Reichardt et al. 2009a, b; Sharp et al. 
2009; Dawson et al. 2006; Sievers et al. 2009). Our ap- 
proach is to find the sources in the ACT maps, mask 
them, and assess the residual contribution with models. 

We find the sources using a matched filter with noise 
weighting derived from the statistical properties of the 
map (Teginark & de Oliveira-Cost a 1998). We identify 
as sources all pixels detected at > 5 a having at least 
3 neighboring pixels detected at > 3a. A 5o* detection 
corresponds to roughly 20mJy in the filtered map. The 
selection criteria are tested on simulated sky maps to 
determine the sample’s purity and completeness. The 
sample is approximately 85% complete at 20mJy and 
~ 100% complete at 50mJv. For fluxes greater than 
20 mJy, the detections have a purity of ~ 95% (that is, in 
a hundred detections approximately five are false). Sim- 
ulations show that at 20 mJy deboosting (Condon 1974) 
is a ~ 9% effect, and at 40 mJy a < 1% effect. 

In the 228 deg 2 area, we detect 108 sources. Of 
these, 105 can be identified with sources in the 
PMN (Wright et al. 1994), SUMSS (Mauch et al. 2003), 
and/or AT20G (Murphy et al. 2009) radio source cata- 
logs. One source, not in these catalogs, is in the 2MASS 
catalog (and thus has IR emission). Two have no previ- 
ously measured counterparts. In a 1' search radius, 17 
sources have both radio and IR. identifications. With the 
exception of the single 2MASS source, all are treated as 
radio sources. 

We fit the flux distribution of the radio sources to three 
models. The distribution agrees with the Toffolatti et aL 
(1998) radio model scaled by 0.49 ±0.1 2 in source counts, 
where the uncertainty is statistical only. In the WMAP 
analysis of sources at 41 GHz, Hinshaw et al. (2007) 
found a good fit with the same model after scaling the 
counts by 0.64. Thus it appears that on average the ra- 
dio population at 150 GHz is composed predominantly 
of flat-spectrum sources. We prefer to compare the data 
to models by scaling the model’s flux rather than the 
counts, as the extrapolation of the radio flux into the 
150 GHz range is typically the most uncertain part of 
any model. For the Toffolatti et al. (1998) model, the 
required flux scaling is 0.5. The sources are also well fit 
by the De Zotti et al. (2005) model scaled by 1.2 in flux. 
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Lastly, we fit the radio sources model in Sehgai et ah 
(2010) and find that its flux should be scaled by approx- 
imately 2.1. All scale factors have a statistical uncer- 
tainty of about 25%. Each of these models can be used 
to estimate the residual rrris flux below the 20mJy cut, 
as discussed in Section 7.1. We plan a more thorough 
exploration of the radio sources in an upcoming publica- 
tion. 

Before computing the power spectrum, we mask a 10' 
diameter region around each of the 108 sources. We call 
this Mask-ACT. In total, 2.3 deg 2 are masked, or 1.0% 
of the map. Doubling the size of the masked holes has a 
negligible effect on the power spectrum. We also create 
a cluster mask, Mask-C, by finding all 5a clusters and 
check that using this mask combined with Mask-ACT 
has a negligible effect on the power spectrum. We also 
generate a mask, Mask-R, based on known radio sources 
measured in the SUMSS (0.8 GHz), PMN (5 GHz) and 
ATCA (20 GHz) catalogs. A model that has been fit to 
multiple data sets (Sehgai et al 2010) is used to estimate 
the mean spectral index, and this index is then used to 
determine the flux cut in the radio catalogs. Since the 
model is based on ensemble properties of radio sources it 
may miss lower-flux low frequency sources having shallow 
indices. The resulting power spectrum is independent of 
which radio source mask we use, or if we use the union 
of the two, and does not depend on whether we include 
the cluster mask. 

Diffuse dust emission becomes dominant near the 
galactic plane. We check for evidence of dust in the 
region of our map nearest to the plane, from right as- 
cension 6 h to 7 h . Through a cross-correlation with 
the estimated dust map of Finkbeiner et al. (1999) (the 
FDS map), we can limit the dust contribution here to 
£(£ + l)Ci/(2n) < 5 pK 2 at £ = 1000, while a direct 
power spectrum of the FDS map in the same region gives 
^ 1 pK 2 . As the diffuse dust component decreases with 
increasing £, it is not significant in our analysis. 

5. POWER SPECTRUM METHOD 

We estimate the CMB power spectrum using the adap- 
tive multi-taper method (AMTM) with prewhitening, 
described in Das et al. (2009). We make independent 
maps from subsets of the data and use only cross-spectra 
between maps to estimate the final power spectrum. All 
operations are performed using the flat-sky approxima- 
tion. We summarize the method in this section. 

5.1. Fields used for power spectrum analysis 

We find the power spectrum of our map by separate 
analysis of each of the 13 patches shown in Figure 2. 
Each patch is 4? 2 x 4? 2 in size, and together they cover 
a rectangular area of the map from a = 0 h 48 m to 6 h 52 m 
(12° to 103°) in right ascension and from S = —55° 16' 
to — 51°05 / in declination. This area is the region of a 
larger survey having the lowest noise. We also split the 
raw data into four subsets of roughly equal size, with 
the data distributed so that any four successive nights 
go into different subsets. The four independent maps 
generated from these subsets cover the same area and 
have approximately the same depth. All maps are fully 
cross-linked. That is, they all contain data taken with 
the sky both rising and setting. 


5.2. Spectrum of a single patch 

We estimate the spectra of the 13 patches indepen- 
dently, before taking a weighted average to find the final 
spectrum. There are four independent sky maps, from 
which six cross-spectra and four auto-spectra are evalu- 
ated on each patch. We use a weighted mean of only the 
cross-spectra for the final spectral estimate. The weights 
depend on both the cross- and auto-spectra, as discussed 
below. 

Before separating the four maps into 13 patches, each 
beam-convolved map, 27(0), is initially filtered in Fourier 
space with a high-pass function F c (£). This filter sup- 
presses modes at large scales that are largely uncon- 
strained. These modes arise from a combination of in- 
strument properties, scan strategy, and atmospheric con- 
tamination. We choose a squared sine filter, given in 
Fourier space by the smooth function 

f 0 

F c (£) - { sin 2 x(£) 


l «min 

Gain F £ <. £ m a,x (2) 

£<£ 


where x(£) ~ (ir/2)(£ ~ 4iin)/(4nax - £ min ). We choose 
£min — 100 and t max = 500. The 13 patches of the map 
are treated separately from this point. The map of each 
patch is prewhitened using a local, real-space operation 
to reduce the dynamic range of its Fourier components 
(Das et al. 2009). The prewhitening operation involves 
adding a fraction (2%) of the map to an approximation 
of its Laplacian. The Laplacian is computed in real space 
by taking the difference between the map convolved with 
disks of radius V and 3k Then the maps are multiplied 
by the point source mask. The prewhitening step greatly 
reduces the leakage of power from low to high multipoles 
caused by the action of the point source mask on the 
highly colored CMB power spectrum. 

For each patch, w r e compute the six 2D cross- spectra 
and four 2D auto-spectra. The axes correspond to right 
ascension and declination. Each spectrum is computed 
using the adaptively weighted multi-taper method, using 
JVtap = 5 2 tapers having resolution parameter N re s = 3 
(see Das et al. 2009) . Windowing the maps with 25 or- 
thogonal taper functions allows us to extract most of the 
statistical power available in the maps, at the expense of 
broadening the resolution in angular frequency by a fac- 
tor of approximately N res . At this stage, each 2D cross- 
spectrum, ^ , between submaps a and fi on patch i 

incorporates the effects of the filter, prewhitening, ta- 
pering, the point source mask, and the beam; they are 
analogous to a “pseudo power spectrum” of an apodized 
map (e.g., Hivon et al. 2002). Each 2D spectrum is then 
averaged in annuli with a narrow range of \t\ to give the 
binned pseudo-spectrum with 

B)f = P bf &f. (3) 

Throughout this paper we define B$ = £.(£-\-l)Ci/2r. The 
binning function Pm is set to one for pixels lying in an 
annular bin (indexed by b) of width A£ = 300 centered on 
\J£\ ~ £b, and zero elsewhere. The size of the patches and 
the resolution of the tapers dictates the width of the bins. 
For our square patches of side s = 4? 2, the fundamental 
frequency resolution in Fourier space is S£ = 2tt/s w 90. 
The application of the tapers with a resolution param- 
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eter jV res = 3 further degrades the resolution to N res 5£, 
so bins chosen to be smaller than At ~ 270 would be 
unavoidably correlated. The binning function is also set 
to zero where \£ x \ < 270. This region of Fourier space is 
particularly sensitive to scan-synchronous effects, either 
fixed to the ground or in phase with the azimuth scan. 

The binned pseudo spectrum B l b a/3 is then deconvolved 
with the mode- mode coupling matrix, which takes into 
account the combined effects of tapering and masking 
and can be computed exactly. Lastly, we divide by the 
£-space representations of the prewhitening filter, of the 
high- pass filter F c {£b ), and of the beam to obtain an un- 
biased estimate of the true underlying spectrum B 1 ^ . 

This procedure has been tested with simulations. The 
number and resolution of the tapers are chosen as the op- 
timal balance between maximizing information and min- 
imizing bias caused by leakage of power. The simulations 
confirm that increasing the number of tapers beyond 5 2 
has a negligible effect on the spectrum errors. 


5 . 3 . Combining patches 

The final power spectrum estimator is given by a 
weighted mean over N — 13 patches. 


B b = 


AU wjBb 

v—V-V j ' 

£<=i < 


(4) 


where 4 B\ = Yla p -,a<p B l b a ^/6 * s ^ le Inean of the six 
deconvolved cross-power spectra in patch i (assuming 
equal weights), arid a and ,8 index the four independent 
maps of that patch. The weights are chosen as the in- 
verse of the variance of this estimator in each patch, i.e. 
wl = l/a 2 (B l h ) ) where 


v 2 (Bi) = (c Blf) - ( Bt ) 2 (5) 


and the average is taken over the several cross-spectra 
computed for patch i. The first term contains the 4-point 
function of the temperature field and is approximated as, 

X (BfBf + Bf B^ + ), (6) 

a,/?, 7 , <5 
a</?;7<<5 


where M — 4 is the number of submaps per patch. 
We have neglected any connected (non-Gaussian) part 
of the 4-point function due to components such as point 
sources. This is a reasonable approximation when choos- 
ing the 13 weights, because the expression (Equation 6) 
is dominated by the auto-spectrum terms, which in turn 
are noise-dominated. 


5 , 4 . Power spectrum covariance 

We estimate the band power covariance matrix E using 
the scatter in the power spectrum among the N ~ 13 
patches, 

£w' = <A B b AB b <) (7) 

4 Here we introduce the notation X as an unbiased estimator of 
the quantity X, in the sense (x'j — X. 



FlG. 3. — The total estimated error a on the power spectrum 
(blue) given by the analytic expression (Equations 9 and 10). The 
uncertainty on a is found using simulations and is shown by the 
shaded blue band. The Gaussian sample variance (red line, labeled 
“S”) dominates for l < 1200, and atmospheric plus instrument 
noise (purple, “N”) dominates at £ > 2000. The non-Gaussian term 
due to unmasked point sources and clusters of galaxies contributes 
about 15% of the variance at 2500 < i < 6000 (green, “P”). The 
errors estimated using the scatter of the results from the thirteen 
patches (Equation 8 are shown for comparison (black points); they 
agree well with the analytic errors. 


E"i w\K,((B\-B b )(%-B v )) 

~ Y~>iV i j ' (®) 

£i=i w l Ei=i K 

The square roots of the diagonal elements of the covari- 
ance matrix are reported as the errors on our power spec- 
trum estimate. 

To test the accuracy of the error estimate, the errors 
on the power spectrum are also computed analytically. 
Three terms contribute: sample variance in the CMB 
multipoles due to limited sky coverage, instrumental and 
atmospheric noise, and a non-Gaussian term due to un- 
masked point sources and galaxy clusters. The diagonal 
component of the variance in one patch can be written 
as the sum of these terms, in order: 


<r 2 (C b ) 


2Cf, 

rib 


A + 4 CbNb/M + N b 2 /n u 


n b 


°p 

/sky 


(9) 


where N b ~ (C& 4 — C&), the difference between the 
weighted mean auto- and cross-spectra, estimates the av- 
erage power spectrum of the noise; M — 4 is the num- 
ber of maps with independent noise properties per patch; 
n w — M (M — 1)/2 = 6 is the number of cross- spectra per 
patch; rib counts the number of Fourier modes measured 
in bin b (that is, the number of pixels falling in the appro- 
priate annulus of Fourier space); f S k y is the patch area 
divided by the full- sky solid angle, 47rsteradians; and C b 
is the weighted mean cross-spectrum. In the last term, 
Op is given by the non-Gaussian part of the four-point 
function. Thus 

4 = T [(ar*(ft)> - 3 ( 5T 2 (h )) 2 ] , (io) 

where T is the temperature map and the average is over 
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pixels in the map. The term arises from the Poisson- 
distributed components in the maps, including unre- 
solved point sources and clusters of galaxies. For pur- 
poses of the covariance calculation, we assume that such 
sources are distributed independently of one other. This 
term is constant as a function of £ and is computed from 
the masked maps with the high-pass filter (Equation 2) 
applied. The overall variance is given by the weighted 
mean of the patch variances (Equation 9). The covari- 
ance among bins is small in the limit that the four-point 
term does not dominate, and with bins chosen to be suf- 
ficiently large. 

Figure 3 shows the three components of the errors. The 
Gaussian part due to sample variance dominates below 
t — 2000; it scales as Cb/ yfnb- The Poisson term due to 
the non- Gaussian clusters and unresolved point sources 
contributes about 15% of the error budget for 2500 < 
t < 6000 but is sub-dominant at all scales. The noise 
term dominates at £ > 2000; atmospheric noise is the 
main contribution at angular scales of £ < 1000. As a 
cross-check of the error estimates, we also plot the errors 
derived from the scatter of the power spectra from the 
thirteen patches. 

The uncertainty in the analytic errors is shown as the 
shaded band in Figure 3. It was estimated by Monte 
Carlo simulations. One thousand patches were simu- 
lated with white noise and their power spectra taken by 
the same methods used on the ACT maps. The results 
demonstrate that the errors estimated from the scatter 
among our 13 patches are consistent with the expected 
uncertainty. The same simulations were used to verify 
that the covariance between different power spectrum 
bins (Equation 8) is less than 1%. 

6. POWER SPECTRUM RESULTS 

The binned estimate of the power spectrum Bb is shown 
in Figures 4 and 1, and bandpowers are given in Table 1. 
With our method the bandpowers are estimated to have 
less than 1% correlation with neighboring bins, but the 
window functions have a small overlap, which we account 
for in the analysis. At £ < 2500, the estimated power 
is consistent with previous observations by ground and 
balloon- based experiments. The features of the acoustic 
peaks are not distinguished with the coarse binning, but 
with the fluctuation band-power measured to 5%, this 
spectrum offers a powerful probe of cosmological fluctu- 
ations at small scales. A clear excess of power is seen at 
£ > 2500 which can be attributed to point sources and, 
to a lesser extent, to the SZ effect. 

The large-scale modes are recovered only after iterat- 
ing the maps. To test for convergence of noiseless maps, 
we compute the spectra Cl of difference maps, between 
the processed simulation map and input simulated sky, at 
successive iteration numbers i. To account for differences 
that are non-uniform between data subsets, we use the 
auto-spectra. We show that the amplitude of fluctuations 

in the difference map (y Cf } ) is small, less than 1% of the 

amplitude of fluctuations in the input map (\fc~h) at all 
scales by iteration 500, as shown in Figure 5. To test 
for convergence in the data, where we do not know the 
input map, we estimate the maximum change in power 
between the processed map at iteration i and the final it- 


eration, estimated as 2 \JCyCb using auto-spectra. Here 

C'l is the spectrum of the difference map between itera- 
tion i and the final iteration, number 1000. We define the 
convergence ratio r c as this change in power given as a 
fraction of the uncertainty in the power, a((%), and find 
it to be sufficiently small (less than 0.5) by iteration 500 
at all scales. The cross-correlation with WMAP suggests 
that the maps are well converged down to £ ~ 200. We 
do not divide the angular power spectrum by a transfer 
function at any value of £. 

We test the isotropy of the power spectrum by estimat- 
ing the power as a function of phase 6 — t&n"~ 1 (£ y /£ x ). 
We compute the inverse- noise- weighted two-dimensional 
pseudo spectrum coadded over map subsets and patches. 
The mean cross-power pseudo spectrum is shown in Fig- 
ure 6, indicating the region masked at \£ x \ < 270. The 
spectrum is symmetric for £ to as it is for any real- 
valued maps. To quantify any anisotropy, the power av- 
eraged over all multipoles in the range 500 < £ < 8000 
is computed in wedges of A# = 20°. It is found to be 
consistent with an isotropic 2D spectrum. 


TABLE 1 

Anisotropy power in bands of width A£ — 300. 
Temperatures are in CMB units. 


Central 

h 

l(t + l)C(,/2ir 

( MK 2)ab 

ACDM 

model c 

SZ 

model d 

Source 

model 6 

Total 

750 

2317 ± 98 

2106.9 

3.6 

0.8 

2111.2 

1050 

1313 ±42 

1204.2 

4.6 

1.5 

1210.3 

1350 

882 ± 20 

777.3 

5.3 

2.4 

785.0 

1650 

450 ± 22 

431.9 

5.8 

3.6 

441.3 

1950 

284 ± 14 

248.3 

6.1 

5.0 

259.5 

2250 

153 ± 12 

133.8 

6.3 

6.7 

146.8 

2550 

77 ±9 j 

72.3 

6.5 

8.6 

87.4 

2850 

65 ±9 

38.0 

6.6 

10.8 

55.3 

3150 

43 ±8 1 

19.8 

6.7 

13.1 

39.6 

3450 

29 ±7 

10.3 

6.7 

15.8 

32.7 

3750 

50 ±7 ! 

5.3 

6.7 

18.6 

30.6 

4050 

47 ±12 

2.9 

6.6 

21.7 

31.2 

4350 

35 ± 14 

1.6 

6.6 

25.0 

33.2 

4650 

34 ± 19 

0.9 

6.5 

28.6 

36.0 

4950 

27 ±22 

0.6 

6.4 

32.4 

39.4 

5250 

77 ±33 

0.4 

6.3 

36.5 

43.1 

5550 

38 ±25 

0.3 

6.2 

40.7 

47.3 

5850 

42 ±30 

0.2 

6.1 

45.3 

51.6 

6150 

26 ±33 

0.2 

6.0 

50.0 

56.2 

6450 

103 ± 41 

0.1 

5.9 

55.0 

61.1 

6750 

-30 ± 68 

0.1 

5.8 

60.3 

66.2 

7050 

209 ± 115 

0.1 

5.7 

65.7 

71.6 

7350 

60 ± 139 

0.1 

5.6 

71.4 

77.1 

7650 

-38 ± 213 

0.1 

5.5 

77.4 

83.0 

7950 

-63 ± 349 

0.1 

5.4 

83.6 

89.1 


a There is negligible covariance between bins. The maximum- 
likelihood fit agrees best with ACDM if the Q, data given here 
are multiplied by 0.91. 

b For comparison with the SPT results of Lueker et al. (2009), 
we also compute the spectrum with wider bins. For five bins of 
width A£ — 400 from 3000 < £ < 5000, we find 39 ± 6, 41 ± 7, 
49 ± 10. 36 ± 12, and 35 ± 13. For four bins of width A£ = 900 
from 5000 < £ < 8600, we find 56 ± 19, 45 ± 25, 140 ± 70, and 
80 ±250 pK 2 . The ACT data are also shown in Figure 1 with bins 
of width A£ — 600 starting at £ = 4200 — 4800; the bandpowers are 
given by 31 ±11, 39±18, 36+19, 61±26, 29±58, and 28±116pK 2 . 

c This and all “model” columns are in the same units as the data: 
thermodynamic pK 2 . 

d Assumes the best-fit value of Asz = 0.63. 

e Assumes the best-fit value of A p = 11.9 and uncorrelated Poisson 
sources. 
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Fig. 4. — The observed power spectrum in bandpowers at 150 GHz from ACT observations (points with error bars). At large angular 
scales there is good agreement with the lensed ACDM model of the primary CMB (light blue curve shown for t < 4700). The x 2 of the 
model is 7.1 for 4 ACT data points in the range 600 < £ < 1800. The best-fitting model to the full dataset is shown (dark blue, the highest 
curve for all t > 2500). The complete model includes the primary CMB model plus both a Poisson power from point sources and SZ power 
from clusters; both additional components have been allowed to vary. The complete model has been smoothed by convolution with a boxcar 
window function of width A£ = 300; the primary CMB model has not been smoothed. The narrower, gold band shows the marginalized 
95% CL limits on the Poisson amplitude, while the curve indicates the best-fit amplitude A p = 11.9 pK . The wider pink band shows the 
95% CL upper bound on the SZ amplitude, Asz < 1-63; the dark curve inside it shows the best-fit value of Agz = 0.63. The Poisson and 
SZ power are consistent with higher frequency observations and with ACDM predictions. The fitting procedure is described in Section 7. 


We test that the signals in separate data subsets are 
consistent by taking the cross-spectrum of two difference 
maps formed from the temperature maps, T% of the four 
data subsets via: 

T 12 (h) = {T 1 (n) - T 2 (n)]/2 (11) 

T 34 (n) = [T z {n) - T 4 (n)]/2. (12) 

(The data subsets are described in Section 5.1.) The 
difference maps are expected to contain noise but no 
residual signal. We estimate the cross-spectrum of the 

difference maps, Cb =• (t 12 T 34 ^ using the methods de- 
scribed in Section 5 for M = 2 data segments. The two 
other permutations of the data, C& — (t is T^ 4 \ and 

/t 14 T t23 \ are also tested. The three difference spectra, 
shown in Figure 7, are consistent with no signal. 

7. CONSTRAINTS ON SZ AND IR EMISSION 

We perform a simple analysis of the power spectrum 
to quantify the combined contribution from dusty galax- 
ies and radio sources, and the level of SZ emission. 
We assume a ACDM cosmology with lensing of the 
CMB and parameters from the five-year WMAP anal- 


ysis combined with BAO and supernovae measurements 
(Komatsu et al. 2009). We defer a full investigation of 
cosmological parameter constraints until we improve the 
absolute calibration and better account for astrophysical 
foregrounds. 

Our model for the power at 148 GHz is 

Bf = Sf MB + AszBf + Ap (^) + Br\ (13) 

where is the lensed primary CMB power spectrum; 
Bf z is a template spectrum corresponding to a prediction 
for the SZ emission in a model with o$ = 0.8 at 148 GHz; 
A p quantifies the Poisson point source power, required 
to be positive; and Bf yrr corresponds to correlated point 
source power from clustered galaxies. The SZ template 
we use includes the correlated thermal and kinetic SZ 
(kSZ) effect derived from numerical simulations and is 
described in detail iri Sehgal et al. (2010). Its amplitude 
is assumed to scale with os as the seventh power, such 
that Asz — (crf z /0.8) 7 for fixed baryon density. This ac- 
counts approximately for the frequency-dependent com- 
bination of the thermal SZ component scaling as the 7.5 
to 8th power of og and the sub-dominant kSZ scaling as 
the fifth power. The expected point source power C p is 
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FlG. 5. — Convergence of the maps as a function of iteration. 
The mapmaking algorithm converges well by iteration 500. For 
signal-only simulations (top), the amplitude of fluctuations in the 
difference between the processed output map at iteration i and the 
input map (denoted (CJ) 1//2 ), is less than 1% of the amplitude of 
fluctuations in the input map (\fCb) at all scales by i = 500. Itera- 
tions i = 5, 20, 100, and 500 are shown. For simulations with noise 
and for the data ( middle and bottom , respectively), convergence is 
tested by estimating the maximum change in power between the 
processed map at iteration i and iteration 1000, as a fraction of the 
uncertainty in the power in the final map. For iteration 500 this 
fraction r c (described in Section 6) is sufficiently small, less than 
0.5 at all scales. 



4 


FlG. 6. The estimated two-dimensional power spectrum Cg. 

multiplied by a factor of t to emphasize the angular, rather than 
the radial, variation. The power is consistent with being isotropic, 
when divided into wedges of A0 — 20°. The vertical lines indi- 
cate the narrow region \£ x \ < 270, where excess power from scan- 
synchronous signals contaminates the power spectrum. This re- 
gion is not used for the power spectrum analysis. The regions near 
(|4 b|, |4/|) « (4000,4000) are more noisy but not biased. A power 
spectrum computed without these regions is consistent with the 
one we present, 

given by 

f'Scut n AT 

Cp = g{vf j S 2 ~ dS (m-K 2 sr), (14) 

an integral over all sources with flux S up to some maxi- 
mum flux S clJ t ; where g(v) — (c 2 /2kv 2 ) x [(e x — l) 2 /x 2 e x ) 



FlG. 7. Cross-spectra of difference maps formed from four 
data subsets. In the top panel, two maps T 12 = T 1 — T 2 and 
T 34 = T s — T 4 are formed; both are expected not to contain signal. 
The cross-spectrum is consistent with no signal. The lower two 
panels show the same cross-spectra for the other two permutations 
of the four data subsets and are also consistent with no signal (x 2 is 
25.2, 27.1, and 28.4 in the three panels with 25 degrees of freedom). 

converts flux density to thermodynamic temperature, 
and x ~ hv/kTcMB- Then A p is the binned £(£-\~1)C p /2tt 
Poisson power at pivot £q = 3000. In thermodynamic 
pK 2 , C p = 0.698 x 10“ 6 A p . The conversion to Jy 2 sr _1 
at 148GHz is Gp^sr" 1 ] - 1.55C*[10“ 5 pK 2 sr]. 

The infrared sources are expected to be clustered. 
At small angles galaxies cluster with typical correlation 
function C(Q) oc 9~~ 0 - 8 (e.g., Peebles 1980), which would 
give C# oc £~ 1 - 2 ori non-linear scales. Motivated by this, 
we first adopt a simple template for the correlated power, 

8 ” ,r - 4 “ (sffi) • < 15 > 

and fit for the amplitude A c * Note that on larger 
scales, t < 300, #| orr is expected to flatten and gradu- 
ally turn over (e.g., Scott & White 1999), but at these 
scales at 148 GHz the CMB dominates. Models sug- 
gest that the power from source clustering is less than 
the Poisson component at scales smaller than t = 2000 
(Scott h White 1999; Negrello et al. 2007; Righi et ah 
2008). This is consistent with observations at 600 GHz by 
BLAST (Devlin et al. 2009; Viero et al. 2009). We there- 
fore impose a prior that the correlated power be less than 
the Poisson power at £ — 3000, i.e., 0 < A c < A p . This 
model is likely too simplistic, and the correlated power 
may have an alternative shape (e.g., Sehgal et al. 2010, 
or the halo model considered in Viero et al. 2009). If the 
power from dusty galaxies is instead better described by 
a linear matter power spectrum scaled by a bias factor, 
B| orr becomes almost degenerate with the SZ component 
Bf (see, e.g.. Hall et al. 2009). Since it is not possi- 
ble to separate these components with data at a single 
frequency, our limit on the SZ contribution at 148 GHz 
should be considered as the upper limit on the sum of 
the SZ power and a degenerate correlated source power. 
Multi-frequency data will enable us to investigate the 
shape and amplitude of the correlated component more 
fully. 

Given this model, the likelihood of the data is given by 
- 2 InJf = (Bf - B b ) T ^ l (Bl h - B b ) + Indet £, (16) 
with covariance matrix S defined in Equation 8. The 
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Fig. 8. — Probability distributions for the amplitude A v of an IR Poisson point source contribution of the form Bg ~ £{t + \)Cg/2'x ~ 
A p (£/ 3000) 2 , and of Asz ♦ the amplitude of the SZ effect relative to that of a model with ag = 0.80. Each panel shows two models for the 
point source power spectrum. In one case, we assume the point sources are uncorrelated Poisson-distributed sources. In the second case, 
we marginalize over the amplitude of a correlated term scaling as Bg <x £, Left: The two-dimensional distribution of A v and Asz* The 
filled (blue) regions assume uncorrelated sources; the unfilled (red) regions allow for correlated sources. Center: The distribution of A p 
marginalized over Asz- I n the center and right panels, the curve peaking at higher amplitude (blue) assumes uncorrelated sources. The 
point source power is consistent with SCUBA and BLAST data at higher frequencies. Right: The distribution of Asz marginalized over 
Ap. The shaded region shows the 95% CL limits on erg inferred from WMAP data combined with distance priors. 


theoretical spectrum Bf 1 is computed from the model 
using &l h — where wte is an approximate form 

of the band power window function in band b. We 
marginalize over the calibration uncertainty analytically 
(Ganga et ah 1996; Bridle et ah 2002). 

The uncertainty on the shape of the window function is 
small, of order 1.5% for the 148 GHz band, compared to 
12% overall calibration uncertainty in power. The win- 
dow uncertainty is therefore neglected throughout this 
analysis. We verify that this approximation is valid by in- 
cluding the window function uncertainty in the likelihood 
calculation. The beam Legendre transform is first ex- 
panded in orthogonal basis functions, and the uncertain- 
ties on the basis function coefficients are used to derive 
the window function covariance matrix. The covariance 
is dominated by a small number of modes, so a singu- 
lar value decomposition is taken. The ten largest modes 
are included in the likelihood, following the method de- 
scribed in Appendix A of Hinshaw et al. (2007). We find 
that including the window" function uncertainty has only 
a negligible effect on parameter estimates (< 0.05a). The 
ACT/MB AC beam measurements and the orthogonal 
function expansion are described in detail in Hincks et al. 
(2009). 

We confirm that the expected values of parameters 
Asz, A p , and A c are recovered from maps of the ACT 
simulated sky. We tested simulations with and with- 
out noise in the map, and simulations with realistic 
timestream noise run through the mapper. 

In the range 600 < £ < 1800, where foreground emis- 
sion and secondary effects are sub-dominant, the data are 
consistent with the lensed A COM model alone, with y 2 
— 7.1 for four degrees of freedom (dof). If we rescale the 
maps to check consistency with the model (multiplying 
temperatures by 0.96, a 0.7a change in the calibration), 
then the y 2 value becomes 3.0 with 3 dof. We do not 
rescale the maps in the analysis that follows, though the 
scale factor is allowed to vary when we marginalized over 
uncertainties. 


Using the full range 600 < t < 8100, we find marginal- 
ized constraints on A p = 11.2 ± 3.3 pK 2 (thus, C p = 
(0.78 ±0.23) xl0~~ 5 pK 2 ) and A SZ < 1.63 (95% CL). 
The minimum y 2 = 27.0 for 23 dof. Assuming the scal- 
ing of Asz as (af^) 7 , this implies an upper limit of 
a| z < 0.86 (95% CL). The one- and two-dimensional 
distributions are shown in Figure 8, with limits given 
in Table 2. Marginalizing over the possible SZ power, 
ACT detects a residual point source component at 3a 
(<5y 2 = 10). Marginalizing over a correlated term with 
A c < A p gives A p ~ 9.7 ± 2.8 pK 2 , and af z < 0.84 with 
y 2 = 26.7 for 22 dof. Because A c is forced to be pos- 
itive, its inclusion has the effect of lowering the limits 
on A p and Asz- Nevertheless, we do not find evidence 
for a correlated component with current sensitivity lev- 
els. The estimated parameters vary by less than 0.6a 
when the minimum angular scale is varied in the range 
5000 < Cnax < 8000, or the SZ template is replaced by 
the spectrum of Komatsu Sz Seljak (2002), which is ap- 
proximately 15% lower than our template in the relevant 
range of £ from 1000-5000. We note that there are mul- 
tiple models for predicting Asz (e.g., Bond et al. 2005; 
Kravtsov et al. 2005) and that the relation between Asz 
and ag is an active area of research. 

We also combine the ACT spectrum with the WMAP 
5-year data (Dunkley et al. 2009) to constrain the six- 
parameter ACDM model (defined by the baryon den- 
sity, cold dark matter density, cosmological constant, 
optical depth to reionization, and the amplitude and 
scale dependence of primordial fluctuations at k = 0.002 
Mpc”" 1 ). We model the SZ and point source contribution 
using Equation 13, neglecting a correlated component. 
We find similar results for the point source and SZ am- 
plitude in this extended model, with A p = 11. 5 ±3. 2, and 
Asz < 1-66 (95% CL). The ACDM parameters are con- 
sistent with WMAP alone, with 100ft&/i 2 ~ 2.27 ± 0.06, 
ft c h 2 - 0.111 ± 0.006, ft A - 0.738 ± 0.030. n, = 
0.964 ± 0.014, r - 0.086 ± 0.017, and 10 9 A* = 2.4 ± 0.1. 
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TABLE 2 

Constraints on secondary anisotropies and extra-galactic foreground emission at 150 GHz. 



^sz 

_sz 

a 8 

A P U = 3000) 

ViK 2 

Cp 

10-5 UK 2 

C p 

Jy 2 sr 1 

X 2 /dof 

Poisson point sources 

< 1.63 

< 0.86 

11.2 ±3.3 

0.78 ±0.23 

1.22 ±0.36 

27.0/23 

Poisson + correlated point sources 

< 1.36 

< 0.84 

9.7 ±2.8 

0.68 ±0.20 

1.05 ±0.30 

26.7/22 


7.1. Comparison to other point source observations 

The residual source level, C p = (0.78 ± 0.23) 

xlO _ 5 fiK 2 , combines power from radio sources and 
dusty galaxies that were not removed by the mask. 
For the three radio source models discussed in Sec- 
tion 4, we expect residual power of A p = 6.4 pK 2 
(TofFolatti et al. 1998, after rescaling by 0.4), A p — 
4.1 pK 2 (De Zotti et al. 2005), and A p = 7 pK 2 

(Sehgal et al. 2010). These models’ amplitudes corre- 
spond to C* adio = 0.43, 0.29, and 0.49 x 10~ 5 pK 2 , re- 
spectively. The correction for those few sources below 
20mJy that we mask can be neglected; it is smaller 
than the spread among the models. Given the uncer- 
tainties in the models, we subtract the typical model 
from the total source level and infer that the component 
from residual IR sources lies in the approximate range 
0.2 < Cf x 10 5 < 1 pK 2 . 

The ACT result on total point source power is similar 
to those of APEX at 150 GHz, which finds C p — Flj^f 
xl(U° pK 2 (Reichardt et al. 2009a), of ACBAR (C p = 
2.7±*;£ xlO~~ 5 pK 2 , Reichardt et al. 2009b), and SPT 
(C p = (0.74 ±0.06) x 10“ 5 pK 2 , Hall et al. 2009). This 
last measurement employs a lower flux level for remov- 
ing discrete sources: 6.4 rriJy versus our cut at 20mJy. 
While this difference means that the C p presented will 
contain roughly three times as much power from radio 
sources as the SPT measurement, we nevertheless find 
that the total residual power due to point sources is con- 
sistent between the two results, given the conservative 
assumption that at least one-quarter of the point source 
power observed by ACT is due to dusty galaxies rather 
than to radio sources. 

The IR source models in the literature for 148 GHz 
make a range of predictions. For example, for sources 
less than 20 mJy the Lagache et al. (2004) model gives 
C p = 40 xlO“°pK 2 , the Negrello et al. (2007) model 
gives C p — 1.2 x!0“° pK 2 , and the Sehgal et al. (2010) 
model gives C p = 17.5 xl0“° pK 2 . The best agreement 
comes with the Negrello et al. (2007) model. 

The emission from dusty galaxies in the rest frame can 
be modeled as Sq(v) x v^B u {T), with emissivity index 
0 ~ 1.5 and Planck function B }y (T). The effective in- 
dex a, where S(v) ex u a (Ce oc v 2ot ) , accounts for the 
redshift of the sources, the intrinsic temperature varia- 
tion T\ and the index 0, SCUBA has observed emis- 
sion at 850 pm (353 GHz), where we expect a similar 
population of galaxies to contribute (e.g., Greve et al. 
2004, 2008). Using SCUBA galaxy number counts and 
a model for dA r / d 5, Scott & White (1999) estimate the 
Poisson power to be C p — 730 x 1 0 ~ 5 pK 2 ( 1 90 Jy 2 sr _1 ) 
for S C nt — 50mJy. Combined with residual IR source 


level observed bv ACT, this implies an effective spectral 
index of <2150-350 between 2.6- 3.3. This is consistent 
with emission from dusty star- burst galaxies at high red- 
shift and in line with predictions by White & Majumdar 
(2004) and Negrello et al. (2007). It also agrees with the 
a = 2.6 ± 0.6 index inferred from source fluxes mea- 
sured with MAMBO (1.2 mm) and SCUBA (850 pm) 
(Greve et al. 2004), and with a = 2.3 measured from 
AZTEC (1.1 mm) and SCUBA (Chapin et al. 2009). Ob- 
servations by BLAST at 500 pm (600 GHz) have Poisson 
power (2.7±0.2) x 10 3 Jy 2 sr~ ] (Viero et al. 2009). Com- 
bining this with the ACT data leads to an estimate of 
the effective index 2.7 < <2150-600 JS 5*6? consistent with 
findings by APEX. The consistency with <2150-350 sug- 
gests that similar populations are being probed at these 
frequencies, although BLAST is sensitive to a lower red- 
shift range than ACT. 

7.2. Comparison to other SZ observations 

The ACT constraints on the SZ power indicate an am- 
plitude of fluctuations uf z < 0.86 (95% CL). This re- 
sult is consistent with estimates that combine the pri- 
mordial CMB anisotropy with distance measures, ag = 
0.81 ± 0.03 (Komatsu et al. 2009), and improves on SZ- 
inferred limits at 150 GHz from Bolocam (<r| z < 1.57) 
(Sayers et al. 2009), Boomerang (crf z < 1.14 at 95% 
CL) (Veneziani et al. 2009), and APEX (cr| z < 1.18) 
(Reichardt et* al. 2009a). We do not see evidence for an 
excess of SZ power, in contrast to lower frequency ob- 
servations by CBI at 30 GHz which prefer a value 2.5<x 
higher than the concordance value (<Xg Z = 0.922 ±0.047, 
Sievers et al. 2009). 

The ACT result is also consistent with the recently 
reported crf z = 0.773 ± 0.025 (Asz = 0.42 ± 0.21) from 
the South Pole Telescope (Lueker et al. 2009). The SPT 
team prefers a form for the correlated point sources that 
is covariant with the SZ template in the £ = 3000 range 
(Hall et al. 2009). If such a form is correct, then the 
A$z we report should be interpreted as an upper limit 
on correlated point sources plus the SZ effect. Analysis 
of the ACT’s 218 and 277 GHz data will shed light* on 
possible forms of the correlated component. 

The ACT results on ag are also consistent with several 
recent studies based on the analysis of ROSAT X-ray 
flux-selected clusters. We note that these X-ray clus- 
ter studies themselves are consistent with measures of 
ag from richness or weak-lensing selected cluster sam- 
ples (see references in the three articles cited here). 
Henry et ah (2009) find <7 8 (fim/0.32)°* 30 - 0.86 ± 0.04 
(for Q m . < 0.32) using cluster gas temperatures measured 
with the ASCA satellite. Vikhlinin et al. (2009) obtain 
<7g(O m /0.25) a47 — 0.813 + 0.013 (stat) ±0.024 (sys) using 
temperatures derived from Chandra observations. Simi- 
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larly, Mantz et al. (2009) find that in spatially flat mod- 
els with a constant dark energy equation of state, EOSAT 
X-ray flux-selected clusters yield O m ~ 0.23 ± 0.04, 
cr 8 = 0.82 ± 0.05. 

After this article was completed, WMAP 7-year mea- 
surements of the brightness of Mars arid Uranus were 
released (Weiland et al. 2010). They suggest that 7\ j ~ 
107 K, which is ~ 4% dimmer than the value used here. 
If the new Uranus temperature is adopted, then all ab- 
solute temperatures and source brightnesses in this work 
would be reduced by 4%, and all Q values by 8%. 

8. CONCLUSIONS 

We have presented the first map and the power spec- 
trum of the CMB sky made using data from the Atacama 
Cosmology Telescope at 148 GHz. With this map we can 
compare to WMAP at degree angular scales and mea- 
sure point sources with a resolution of 0?02. With an 
unbiased estimator, we extract the power spectrum, CV, 
over a range of power exceeding 1() 4 . 

We have interpreted the spectrum with a simple model 
composed of the primary CMB, a possible SZ contribu- 
tion, and uncorrelated point sources. This analysis pro- 
vides a new upper bound on the SZ signal from clusters 
(erf z < 0.86 at 95% CL, though this is subject to un- 
certainty in the SZ models). A coordinated program of 
X-ray, optical, infrared, and millimeter- wavelength ob- 
servations of the largest SZ clusters is underway. 

These high angular resolution measurements probe 
the microwave power spectrum out to arerninute scales. 
Above £ ~ 2500, the spectrum is sensitive to non-linear 
processes such as the formation of galaxy clusters and 
dusty galaxies. On the \cw-£ end, the spectrum mea- 
sures the Silk damping tail of the CMB which can be 
computed using linear perturbation theory as applied to 
the primordial plasma. It is clear that to understand the 
£ > 1000 end of the primary CMB, and thus to improve 
significantly on measurements of the scalar spectral index 
and its running, source modeling will be required. Future 
analyses of the ACT data will include the two higher- 
frequency channels and additional sky coverage. In an- 
other approach, one can measure the high-f E- modes, 
because the polarized CMB to foreground ratio is ex- 
pected to be higher than that for the temperature. We 
are pursuing both programs, as a polarization-sensitive 
camera is currently under development. 
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